Dynamic and static density-density correlations in the one-dimensional Bose gas: 

exact results and approximations 
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We discuss approximate formulas for the dynamic structure factor of the one-dimensional Bose 
gas in the Lieb-Liniger model that appear to be applicable over a wide range of the relevant param- 
eters such as the interaction strength, frequency, and wavenumber. The suggested approximations 
are consistent with the exact results known in limiting cases. In particular, we encompass exact 
edge exponents as well as Luttinger liquid and perturbation theoretic results. We further discuss 
derived approximations for the static structure factor and the pair distribution function g(x). The 
approximate expressions show excellent agreement with numerical results based on the algebraic 
Bethe ansatz. 
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I. INTRODUCTION 

Correlations in ultracold atomic gases arise from the 
interplay of quantum statistics, interactions and thermal 
and quantum fluctuations. Recently, a lot of progress 
has been made experimentally to probe and character- 
ize these correlations [l], [U, S 13 ■ The one-dimensional 
Bose gas is a particularly interesting system as quan- 
tum correlations generally play a larger role compared to 
three dimensional Bosc-Einstcin condensates and regimes 
with very different correlation properties can be probed 
experimentally 0, Q- In these experiments, elongated 
"spaghetti" traps are created by optical lattices, which 
confine the atomic motion in the transverse dimensions 
to zero-point quantum oscillations [7|. Thus, the systems 
become effectively one-dimensional. 

Theoretically, interactions of the rarefied atoms in one- 
dimensional waveguides are well described by effective 
S- function interactions Q . The resulting model of a one- 
dimensional Bose gas is an archetype of an integrablc 
but non-trivial many-body system that has been receiv- 
ing long standing interest from physicists and mathemati- 
cians alike. The model was first solved with Bethe ansatz 
by Lieb and Liniger [§, , who calculated the ground- 
state and excitation energies. Depending on the value 
of the dimcnsionless coupling strength, the Lieb-Liniger 
model describes various regimes with the corresponding 
correlations. Being exactly solvable, the model, however, 
does not admit complete analytic solution for the cor- 
relation functions. Up-to-now, this is complicated and 
challenging problem in ID physics [ill, [T3|. 

Dynamical density-density correlations can be mea- 
sured in cold atoms by the two-photon Bragg scattering 
[H, EH • Theoretically, they are described by the dynamic 
structure factor (DSF) [lj] 



S(k,u) = L 



dtdx 
2nh 



,i(ujt—kx) 



(0\Sp(x,t)5p(0,0)\0), (1) 



fluctuations and n = N/L is the equilibrium density of 
particles. We consider zero temperature, where (0| . . . |0) 
denotes the ground-state expectation value. The DSF 
is proportional to the probability of exciting a collective 
mode from the ground state with the transfer of momen- 
tum k and energy Hu>, as can be seen from the energy 
representation of Eq. JTJ) 



S(k,w) = mSp k \n)\ 2 S(ha; - E m + E ) 



(2) 



where 5p{x, t) = p(x, t) — n is the operator of density 



where 5pk = Ylj e j ~ NA(k) is the Fourier compo- 
nent of Sp(x), A(fc) = 1 at k — and A(fc) = otherwise. 
Once the DSF is known, the static structure factor S(k) 
and the pair distribution function g{x) can be calculated 
by integration as is discussed in Sec. III.B. 

Previously known results for the DSF of the one- 
dimensional Bose gas come from Luttinger liquid the- 
ory, which predicts a power-law behavior of the DSF 
at low energies in the vicinity of the momenta k = 
0, 27m. Aim . . . and yields universal values for the expo- 
nents [15|, [H, In the regime of strong interactions, 
we have previously derived perturbatively valid expres- 
sions covering arbitrary energies and momenta at zero 
[IH and finite temperature For finite systems, it 

is possible to compute the correlation functions numer- 
ically, using the results of the algebraic Bethe ansatz 
calculations [13, Hl|. Finally, the exact power-law be- 
havior along the limiting dispersion curve of the collec- 
tive modes has recently been calculated in Refs. [2^.[23j. 
These exponents differ from those predicted by Luttinger 
liquid theory raising the question whether the different 
results are compatible with each other. We address this 
question in Sec. IIII Al of this paper, where we show that 
the results can be reconciled by taking appropriate lim- 
its. The apparent difference between the edge exponents 
valid along the dispersion curves and the Luttinger liquid 
result in the limit of vanishing energy can be traced back 
to the fact that the dispersion relations are curved and 
not straight, as is presumed by Luttinger liquid theory. 
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The exact values of the exponents found in Refs. [221 , 
|23| are of importance; however, they are not sufficient 
for practical estimations of the DSF as long as the pref- 
actors are not known. In this paper we construct an 
approximate formula for the DSF [24j based on the ex- 
ponents of Refs. [13, HH. Within the proposed scheme, 
the prefactor can be found using the well-known /-sum 
rule (see, e.g. [HI.) The result turns out to be consis- 
tent with numerical results by Caux and Calabrese [2(| ■ 
Besides, it is compatible with the results of Luttinger liq- 
uid theory [TEHEES] and perturbation theory pjf. The 
approximate formula, in effect, takes into account single 
quasiparticle-quasihole excitations but neglects multipar- 
ticle excitations. We also present an approximate expres- 
sion for the static structure factor and for the density- 
density correlation function, which is derived from the 
approximation for the DSF. 



II. EXACT RESULTS FOR DYNAMIC 
STRUCTURE FACTOR IN LIEB-LINIGER 
MODEL 

We model cold bosonic atoms in a waveguide- like micro 
trap by a simple ID gas of N bosons with point interac- 
tions of strength <7b > 



H 



h 2 d 2 



N 

^~ 7 2~mdx 2 

l—l 1 



9b 5(xi-Xj) (3) 



and impose periodic boundary conditions on the wave 
functions. The strength of interactions can be measured 
in terms of dimcnsionlcss parameter 7 = rag^/ (h 2 n). In 
the limit of large 7, the model is known as the Tonks- 
Girardeau (TG) gas. In this limit, it can be mapped 
onto an ideal Fermi gas since infinite contact repulsions 
emulate the Pauli principle. In the opposite limit of small 
7, we recover the Bogoliubov model of weakly interacting 
bosons. 



A. DSF expansion in I/7 

For finite 7, the model can also be mapped onto a 
Fermi gas [25| with local interactions, inversely pro- 
portional to gs [H, HE H3- Using the explicit 
form of the interactions, one can develop the time- 
dependent Hartree-Fock scheme [TH, [l!| in the strong- 
coupling regime with small parameter I/7. The scheme 
yields the correct expansion of the DSF up to the first 
order [HGH 



S(k,uj) 



4fc 



1 



1 



In- 



O 



1 



N 4fc \ 7/ 27' w^-W \7* 

n (4) 

for uj-(k) ^ u> ^ uj + (k), and zero elsewhere |28[. The 
symbol 0(x) denotes terms of order x or even smaller. 
Here uj±(k) are the limiting dispersions [29j that bound 



quasiparticle-quasihole excitations (see Fig. [T]); in the 
strong-coupling regime they take the form 

uj±(k) = h\2k F k ± fc 2 |(l - 4/ 7 )/(2m) + 0(l/7 2 ). (5) 

By definition, kp = nn and £p = fi 2 fcp/(2m) are the 
Fermi wave vector and energy of TG gas, respectively. 



B. Link to Luttinger liquid theory 

Luttinger liquid theory describes the behavior of the 
DSF at low energies for arbitrary strength of interactions 
[H, QjJ. In particular, one can show 0, H3] that m the 
vicinity of "umklapp" point (k = 2im, uj = 0) it is given 

by 



S(k,i 



N 



2 A' 



L0 2 _(k) 



K-l 



(0) 



for uj ^ uj- (k), and zero otherwise. Within the Luttinger- 
liquid theory, the dispersion is linear near the umklapp 
point: uj-(k) ~ c\k — 2im\. By definition, 



K — him I (mc) 



(7) 



and c is the sound velocity. For the repulsive bosons, 
the value of parameter K lies between 1 (TG gas) and 
+00 (ideal Bose gas). In the strong-coupling regime, the 
linear behavior of the dispersions ([5]) at small momen- 
tum determines the sound velocity, which allows us to 
calculate the value of the Luttinger parameter 



A' = l + 4/ 7 + 0(l/7 2 ). 



(8) 



The coefficient A(K) is model-dependent; in the Lieb- 
Liniger model, it is known in two limiting cases: A(K) = 
tt/4 at K = 1 and A{K) ~ exp(-2 lc K)TT 2 /T 2 {K) 

for K > 1 [13, where 7 C = 0.5772 ... is the Euler con- 
stant and T(K) is the gamma function. 

By comparing the first-order expansion (f5|) in the vicin- 
ity of the umklapp point with Eq. ^ and using the ex- 
pansion (|8|), one can easily obtain the model-dependent 
coefficient at large but finite interactions when A'-l C 1 



A(K) = —[1 — (1 + 4 In 2) (K — l)} + ((K I) 2 ) 



(9) 



Note that the relation ([6]) leads to different exponents 
precisely at the umklapp point and outside of it: 

w 2 ^" 1 ), fc = 27m, 

(uj-uj„) k -\ k^2im. 



S(k,u) 



(10) 



Exact edge exponents from the Lieb-Liniger 
solutions 



As was shown in Refs. [22j,|23j (see also [30j), within the 
Lieb-Liniger model the DSF exhibits the following power- 
law behavior near the borders of the spectrum uj± (k) 



S(k,uj) ~ \oj — uj±(k) 



=F/x±(fe) 



(11) 
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FIG. 1: (Color online) Numerical values of the DSF ([2} for 
the coupling parameter 7 = 10 [20|. The dimensionless value 
of the rescaled DSF S(k,u})eF/N is shown in shades of gray 
between zero (white) and 1.0 (black). The upper and lower 
solid (blue) lines represent the dispersions u>+(k) and ui-{k), 
respectively, limiting the single "particle-hole" excitations in 
the Lieb-Liniger model at T — 0. The dispersions are ob- 
tained numerically by solving Lieb-Liniger's system of inte- 
gral equations (see Appendix [S} . The gray scale plot of the 
DSF demonstrates that the main contribution to the DSF 
comes from the single particle-hole excitations, lying inside 
the region cj-(k) ^ ui sj ui+(k) (see also Fig. [3j) . 



in 10-k plane, where the exponents are different, vanishes 
in the point k = 27m; hence, the Luttingcr exponent 
2 (if — 1) becomes exact there. 

A typical behavior of the exponents is shown in Fig. [2j 
As described in Appendix YK[ the exponents can be easily 
evaluated by solving equation (|A6|) for the shift function 
and using Eq. (|A"9j) . 

D. Algebraic Bethe ansatz 

Recent progress in the computation of correlation func- 
tions within the Lieb-Liniger model and other ID mod- 
els has been achieved through the algebraic Bethe ansatz 
pp| | . In this method, matrix elements of the density op- 
erator involved in Eq. @ were calculated with the alge- 
braic Bethe ansatz. They are given by the determinant 
of a matrix, which can be evaluated numerically for a 
finite number of particles. So, this method is based on 
combining integrability and numerics. The results of the 
numerical calculations of Ref. [2(| are shown in Figs. [1] 
and [51 

III. APPROXIMATE EXPRESSION FOR 
DYNAMIC STRUCTURE FACTOR 
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FIG. 2: (Color online) Typical behavior of the exact expo- 
nents in Eq. JTTJ . The diagram shows fj,± for 7=10 obtained 
numerically using the method of Ref. [U described in Ap- 
pendix [X] 



The positive exponents fj,± [29| are related to the quasi- 
particle scattering phase and can be calculated in the 
thermodynamic limit by solving a system of integral 
equations [23|. In particular, Imambekov and Glazman 
[231 ] found the following right limit 



lim H-{k) 

k — >27rn~ 



2Vk(Vk-i), 



(12) 



which is different from the Luttingcr liquid exponent 
(fit))) . However, Imambekov's and Glazman's result (fT2)) 
is accurate in the immediate vicinity of w_ provided that 
the finite curvature of a->_ (fc) is taken into consideration. 
Thus the difference in the exponents can be attributed 
[23d to the linear spectrum approximation within the Lut- 
tinger liquid theory. Note, however, that the thin "strip" 



A. Approximate expression for arbitrary values of 
interaction strength 

Here we suggest a phenomenological expression, which 
is consistent with all the above-mentioned results. It 
reads 



S(k,u) = C 



(lu c 



')/*+ 



(13) 



for aj_(fc) ^ lu ^ uj + (k), and zero otherwise. It follows 
from energy and momentum conservation that S(k,cu) is 
exactly equal to zero below LU-(k) for ^ k ^ Inn. In 
the other regions of lu > to + and to < w_ (for k > 27m) , 
possible contributions can arise due to coupling to multi- 
particle excitations [lj|. However, these contributions 
are known to vanish in the Tonks- Girardeau (7 — > 00) 
and Bogoliubov (7 — ► 0) limits and are found to be very 
small numerically for finite interactions (2fjj . 

In Eq. (| 13[) C is a normalization constant, fi+(k) and 
/i_(fc) are the exponents of Eq. (fTTj) , and a = 1+ \j\[K. 
From the definition of A' ([7|), one can see that for re- 
pulsive spinlcss bosons K ^ 1, and, hence, 1 < a ^ 2. 
The normalization constant depends on the momentum 
but not the frequency and can be determined from the 
/-sum rule [l4j 



+00 



dio ioS(k, to) = N 



2m 



(14) 



In Eq. (|13[) we assume that the value of the exponent 
fjL-(k = 2im) coincides with its limiting value (|12[) in 
vicinity of the umklapp point. 
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FIG. 3: (Color online) The Dynamic Structure Factor (DSF) 
in the thermodynamic limit. The proposed approximation 
(|13[) (line) is compared to numerical data from Caux and Cal- 
abrese [201 ] (open dots) . The dashed (red) line shows the data 
of Eq. (|13fl convoluted in frequency with a Gaussian of width 
0.07ep/ft in order to simulate smearing that was used in gen- 
crating the numerical results of Ref. [2C|. The numerical data 
of Ref. [2(| suggest that contributions from multi-particle ex- 
citations for uj > (sharp line in parts a and b) are very 
small. Such contributions are not accounted for by the for- 
mula (I13|l . Insert: DSF at the umklapp point in logarithmic 
scale. The graph shows that the DSF behaves as predicted by 
the Luttinger liquid theory (|16|1 with the exponent 2(K — 1), 
where K = 1.402 ... at 7 = 10. 



The most general way of obtaining tu±(k), n±(k), and 
K is to solve numerically the corresponding integral equa- 
tions (see Appendix [A"]) . 

Note that the sum rule for the isothermal compress- 
ibility [3 



lim 

fc^O 



H-oc 



S(k, uj) da; 
Nu> 



2mc 2 



(15) 



is satisfied by virtue of Eq. (|A10[) and the phonon behav- 
ior of the dispersions at small momentum: uj±(k) ~ ck 
(see Fig. Q]). 

Now one can see from (TT3l) that 



S(k,u) 



,2(K-1) 



27m, 



(lu -cu^-( k \ k 7^ 2nn. 



(16) 



Thus, the suggested formula (|T3[) is consistent with 
both the Luttinger liquid behavior at the umklapp point 
and Imambckov's and Glazman's power-law behavior in 
vicinity of it, as it should be. 

In the strong-coupling regime, Eq. (fl3|) yields the cor- 
rect first order expansion ([?]). In order to show this, it is 
sufficient to use the strong-coupling values of K , the 
exponents (|A12|) . and the frequency dispersions ([5]). 

Comparison with the numerical data by Caux and Cal- 
abrese [2(| (Fig. [3]) shows that the suggested formula 
works well in the regimes of both weak and strong cou- 
pling. 

Let us discuss how the Bogoliubov approximation 
arises in the weak-coupling regime in spite of the ab- 
sence of the Bose-Einstein condensation in one dimension 
even at zero temperature [U HU- At small 7, the up- 
per dispersion curve w+(fc) is described well [Hj by the 
Bogoliubov relation [331 ] 



T 2 + 4T k e F -//ir 2 , 



(17) 



where Tk = h 2 k 2 /(2m) denotes the usual one-particle ki- 
netic energy. Besides, when q is finite and 7 — > 0, the as- 
sociated exponents fx + approach the limiting value (|A11[) , 
which in turn is very close to one. This implies that the 
DSF has a strong singularity near u + , and, hence, it is 
localized almost completely within a small vicinity of the 
upper branch (see Fig. [4]). Thus, the behavior of the 
DSF simulates the (^-function spike. One can simply put 
<S f Bog(fc,w) = C5(u> — LUk) and determine the constant C 
from the /-sum rule (|14[) 

S BoK {k,w) = Np-5(u-u k ). (18) 



B. Simplified analytic approximation for 
intermediate and large strength of interactions 

One can further simplify the expression for the DSF 
and replace the parameter a in Eq. (|13p by its limiting 
value a = 2 for the Tonks- Girardeau gas, which turns 
out to be a good approximation even for intermediate 
coupling strength 7^1. This replacement allows us to 
write down the normalization constant explicitly. From 
the /-sum rule we obtain 



S(k,Lu) =N 



k 2 r(2 + ah 



M _) {^-uiy- 



mr(i + /i-)r(i-// + ) {uj 2 + 



,2 s ) AM- 



x (wi - or.) 



2 \fl + -fl- 



(19) 
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FIG. 4: (Color online) Comparison of the two approximations 
for the DSF. The solid (blue) line represents the "universal" 
approximation (|13[) . The dashed (red) line is the approxima- 
tion (|19[) . The two curves coincide almost everywhere except 
for the umklapp point (u> = 0, k = 2im). The "universal" ap- 
proximation reproduces the correct power-low behavior of the 
Luttinger liquid theory S(k, w) ~ lu 2(k ' 1) with K = 3.425 . . . 
at 7 = 1, see the insert. However, the difference in absolute 
values is negligible due to a strong suppression of the DSF 
outside the close vicinity of the upper brunch. 



for aj_(fc) ^ w ?5 Lu + (k), and zero otherwise. This ap- 
proximation ensures all the properties of the DSF men- 
tioned in Sec. HH except for the Luttinger liquid theory 
predictions in the close vicinity of the umklapp point (see 
discussion in Sec. Ill Bj ). However, outside the umklapp 
point, it agrees well with the Caux and Calabrcse numer- 
ical data (see Figs. [3] and 2]). 

From the explicit formula (fT!?)) one can find analytic 
expressions for the static structure factor and the dy- 
namic polarizability. The static structure factor S(k) = 
(pkp-k)/N contains information about the static corre- 
lations of the system and it is directly related to the pair 
distribution function [14, 01 



l(x) = 1 + 



r+°° au 

/ — cos(fcr) \S(k) - ll . 
Jo ™ 



(20) 



The static structure factor can be obtained by integrating 
the DSF over the frequency 



h r +oa 



S(k, uj) duj. 



(21) 



Note that the "phonon" behavior of both dispersions en- 
sures the correct behavior of the static structure factor 
at small momentum. Indeed, it follows from the gen- 
eral expression (fT3"]) that S(k) ~ Hk/(2mc). In the large- 
momentum limit, we have uj + /lj^ ~ 1, which leads to the 
correct asymptotics S(k) — > 1 as k — ► +oo. Equations 



IH) and J5IJ) yield 



S(k) 



=2-Fi( 2+/ i -~A i +, 1 + M-,2 + /i--^ 



hk 2 



2mui+ \uj + 



UJ _ x 1+2 M - 



(22) 



where 2-F1 is the hypcrgcomctric function. The results 
for the static structure factor are plotted in Fig. \5\ One 
can see that the formula for the static structure function 
works well even for weak coupling. This is due to the 
smallness of the DSF contribution to the static structure 
function at the umklapp point for small 7. Thus, the ap- 
proximate formula provides a good accuracy for arbitrary 
strength of interactions. In the weak-coupling regime, 
one can obtain a good approximation for the static struc- 
ture function from the Bogoliubov formula (TT5)) for the 
DSF 



s(k) = A. 

fuj k 



(23) 



The behavior of the pair distribution function (|2"tJ)) in 
the Lieb-Liniger model was studied at large [ll[ and short 
distances jl9l [35L HH in various regimes. For 7 < 1, 
one can obtain from Eqs. ([20)) and (|23|) the analytical 
expression [36| 

g(x) = 1-Vt[£-i(2V7 W*Wi(VrW*-)]> (24) 

where L-i(x) is the modified Struve function and I\(x) 
is a Bessel function. In the opposite limit 7 > 1, one 
can directly use the strong-coupling expression (j4|) for 
the DSF and obtain fijl 



1 



sin 2 z 



2d_ 

7 dz 



z- 
sin z 



2ir d sin 2 z 

7 dz z 2 



4 sin 2 z 

7 z 2 



z 



1 + T) 

dij sin(?7z) In 

1 - ?/ 



+ o( 7 - 2 ), 



(25) 

where z = k^x = imx. The last equation implies that 
g(x = 0) vanishes not only in the TG limit but also in 
the first order of 7~ 1 , which is consistent with the results 

of Refs. [jilajj. 

The behavior of g(x), obtained from the formula (|22j) 
is shown in Fig. [6] It is consistent with both the weak- 
and the strong-coupling limits. 

The dynamic polarizability determines the linear re- 
sponse of the density to an external field [3, H3 ■ It can 
be calculated using the DSF 



2w'S(k,u'] 
w' 2 - z 2 



X{k,z) 

On substituting Eq. (pU in ([26]). we get 

X(k,z) 



doj' 



(26) 



=2*1(1,1 

x N 



M-,2 
1 



(27) 
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FIG. 5: (Color online) The static structure factor versus 
wavenumber for different values of the coupling constant 7. 
(a) The numerical data by Caux and Calabrese [20] (open cir- 
cles) are compared with the proposed analytical formula (|22|l 
(solid lines). The dashed (red) line shows the static structure 
factor in the Bogoliubov limit (|23p . (b) The static structure 
factor obtained with Eq. (|21[) from the general formula for the 
DSF (|13|) is shown by the solid line. These data are consistent 
with the analytical formula (|22|l (dashed line). This indicates 
that the analytical formula for the static structure factor can 
be used even for small values of 7. 

For a retarded response, we should put here z = cj + is. 
At zero temperature the relation S(k,u) = Imx(fc,w + 
ie)/w holds. 

The obtained relations (f2"2")l and (|2"T|) successfully re- 
produce the Tonks-Girardeau limit, considered in detail 
in Refs. [HG3]. 
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FIG. 6: (Color online) The pair distribution function g(x) ver- 
sus the distance in dimensionless units 1/kv (&f = 7rn) for dif- 
ferent values of the coupling constant 7. The solid (blue) line 
represents the pair distribution function (|20p obtained using 
the approximation (|22[l . The Bogoliubov approximation (|23|l 
is indicated by dashed (red) line, while the strong-coupling 
approximation (|25|) by dotted (green) line. The values of g(Q) 
are consistent with the results of Ref. [35j . 



small, anyway, outside the bounds given by the disper- 
sion curves u±. Our formula is consistent with the pre- 
dictions of the Luttinger liquid theory. It gives the exact 
exponents at the edge of the spectrum, the correct first- 
order expansion in the strong-coupling regime, and shows 
good agreement with the available numerical data. For 
intermediate and large values of the interaction strength 
7 > 1 and outside the close vicinity of the umklapp point 
(u> = 0, k = 27m), the further simplified analytic formu- 
las for the DSF (fl9|) and the dynamic polarizability (|27[) 
provide excellent accuracy. The analytic expression (|2"2"j) 
for the static structure factor works well even for weak 
interactions. 

Our results provide a reference against which experi- 
mental measurements of static and dynamic density cor- 
relations in the one-dimensional Bose gas can be tested. 
They further provide a basis for future work on the con- 
sequences of correlations in this interesting system. 
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APPENDIX A: THE LIEB-LINIGER SOLUTIONS 
AND THE EDGE EXPONENTS 

For the Lieb-Liniger model ([3]) of TV bosons embed- 
ded in the circuit of length L, the Bcthe ansatz yields 
[1, EL H3] the following system of equations for TV quasi- 
momcnta A, 



2im-j 



(Al) 



fe=i 



where 8(X) = 2arctan[fi, 2 A/(<7B"i)]; and rij arc integer 
for odd TV and half-integer for even TV. In the same way 
as for the gas of non-interacting fermions, eigenstates of 
Hamiltonian are completely described by a certain set 
of TV numbers rij] the total energy and momentum of the 
states arc given by 



E 



N 



2m 



N 



(A2) 



£4 P N = H^Xj = 2 ^ 

One can also write down explicitly [H, [ill, H3 rather 
cumbersome expression for the cigenfunctions in terms 
of quasimomenta. One can prove the following proper- 
ties of the Lieb-Liniger solutions [HI, HI] : 

i) The Bethe ansatz solutions form a complete set of 
eigenf unctions. 

ii) There is one-to-one correspondence between rij and 
Xj; if rij > nk then \j > Afc. 

Hi) The set of numbers rij = j — 1 — (TV— 1) /2, 1 ^ j ^ TV 
provides the ground state; for the associated quasimo- 
menta we have |Aj| ^90 *S ftw>- 

iv) In the thermodynamic limit, when TV — > oo and 
n = N/L = const, the distance between the consecutive 
quasimomenta tends to zero: Xj — Afc ~ l/L. 

In order to pass to thermodynamic limit, it is conve- 
nient to introduce the density of quasimomenta p(X) = 
lira l/[X(A,+i — Xi)]. From Eq. (|A1[) and property Hi) 
one can obtain Q a system of linear integral equations 
for the density in the ground state 



/go 
K{v, p)p(p) d/j,, 
-an 



'/(I 



p[y) dv. 



i Hi 



We denote 



2,9b 

h 2 



f9Bm\ 2 s 



(A3) 
(A4) 

(A5) 



and assume that quasimomenta fill "Fermi segment" with 
bounds ±<7o- The density at the bounds can be re- 
lated [ll| to the Luttinger liquid parameter by equation 
p(±q ) = VK/(2tt). 

Classification of all the excitations can be done in the 
same manner as for the Tonks-Girardeau gas of non- 
interacting fermions (<7b - ► +oo). Thus, in order to cre- 
ate elementary particle-like excitation, one needs to add 



a quasimomcntum A p beyond "Fermi segment" [— 9o,9o]- 
By contrast, for a hole-like excitation, one needs to re- 
move a quasimomcntum lying inside the "Fermi seg- 
ment". All the excitations can be constructed from the 
above elementary excitations. Adding or removing a par- 
ticle with quasimomcntum X p leads to a new set of quasi- 
momenta Vj . One can conveniently define a shift function 
FB(vj\Xp) = ±(^j — — Vj), where we take plus 

for a particle-like excitation and minus for a hole-like 
excitation. It follows from (|A1[) that the shift function 
obe ys t he integral equation in the thermodynamic limit 



7T + 6(v — X) 

2^ ' 



obeys tJ 
011 

1 f qa 

Fb(v\X)-—J_ K(u,n)F B (n\\)dn 

(A6) 

Using Eq.l (|A2|> and the definitions of the density and 
shift function, one can evaluate the frequency and the 
wavenumber of elementary single excitation in the ther- 
modynamic limit [ll| 



w Pjh (r/) = ± — 
2m 



fep,h(f) = ± 



v 1 - 



2pF B {p\v) dp, 



'la 



'In 



8{v - p)p(n) dp 



'/ii 



(A7) 
(A8) 



For the particle excitation, as discussed above, we should 
take v > 90, and for the hole excitation \v\ < 90. Note 
that Eq. (|A8j) yields k h (±q Q ^0) = T^f and fc p (±9 ±0) = 
±fcp, as it should be. 

Once equations (|A3[) . (|A4|) . and (|A6|) are solved, we 
can write down the dispersions u±(k) [291 ] bounding sin- 
gle particle-hole excitations (see Fig. [T]). The upper lo + 
branch is constructed by creating hole with quasimomen- 
tum 90 and particle with quasimomcntum A > 90, while 
the lower w_ branch by creating hole inside the "Fermi 
segment" (| A| < 90) and particle with quasimomcntum ly- 
ing just above 90. Explicitly, the branches and associated 
momenta corresponds to the pair excitations A/j = Ao 
and A p = A + Ao (upper, or "+", branch), Xh = Ao — A 
and A p = Ao (lower, or branch). Summing energies 
(|A7[) and wavevectors (|A8[) of the corresponding elemen- 
tary excitations yields the dependencies o>(A) and A: (A) 
for the upper branch 

oj + (X) = uj h (X ) + uj p (X + Ao), 
fc+(A) = fc /l (A ) + fcp(A + A ), 

and the lower branch 

w_ (A) = w/j(A - A) + lo p (Xo), 
fc_(A) = fc fc (Ao- A) + fc p (A ). 

After excluding quasimomcntum, we obtain the disper- 
sions as functions of k. 

It turns out [ll[ that the shift function in the thermo- 
dynamic limit is directly related to the scattering phase 
of two particles (holes) in the presence of filled "Fermi 



8 



segment" . As was shown by Imambekov and Glazman 
[231 ], the edge exponents, specifying the behavior of the 
DSF near the dispersion curves, can be written down in 
terms of the associated phases 5±(A) = 2ttFb{^Qq, A) 

-HU^)'-K^)'-(«) 

The exponents obey [23| the limiting relations 



In the strong-coupling limit, one can show [221 ] that 



H±(k) = 2fc/(7m 7 ) + 0(l/ 7 2 ). (A12) 



/i±(0) = 0, (A10) 
M+(+oo) = 1 - 1/(2K). (All) 
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